function [model, Ypred, blad] = rownanieDzPubPat(Y, U1, U2, na, nb1, nb2, nc, nk1, nk2)

%Y - Patenty
%U1 - medycyna
%U2 - immunobiologia

%U1 - ok
%U2 - ok
Y = [Y(1);Y(1);Y(1);Y];

%kor(1) = KorelacjaPrzes(U, Y, nk);
%if (nb>1)
%    kor(2) = KorelacjaPrzes(U,Y, nk+1);
%end

DAT = iddata(Y, [U1 U2], 1);
model = armax(DAT, 'na', na, 'nb', [nb1 nb2], 'nc', nc, 'nk', [nk1 nk2]);
DATpred = predict(model, DAT, 1);
Ypred = DATpred.y;
blad = SumaKwadratowBl(Y, Ypred);